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Abstract Via the magnification bias, gravitational lensing by large-scale structures causes angular cross-correlations between 
distant quasars and foreground galaxies on angular scales of arc minutes and above. We investigate the three-point cross- 
correlation between quasars and galaxy pairs measurable via the second moment of the galaxy counts around quasars and show 
that it reaches the level of a few per cent on angular scales near one arc minute. Combining two- and three-point correlations, a 
skewness parameter can be defined which is shown to be virtually independent on the shape and normalisation of the dark-matter 
power spectrum. If the galaxy bias is linear and deterministic, the skewness depends on the cosmic matter density parameter 
Oo only; otherwise, it can be used to probe the linearity and stochasticity of the bias. We finally estimate the signal-to-noise 
ratio of a skewness determination and find that around twenty thousand distant quasars e.g. from the Sloan Digital Sky Survey 
should suffice for a direct measurement of Qo- 
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1. Introduction 

It is widely believed that structures and galaxies in the Universe 
formed from gravitational growth of Gaussian primordial mass 
density fluctuations dominated by dark matter. Direct support 
for this picture is provided by the recent weak lensing surveys 
of galaxies which measured the systematic distortion of faint 
background-galaxy images produced by the gravitational tidal 
field of intervening dark-matter inhomogeneities : the cosmic 
shear (Bacon et al. 2000 and 2002; Hammerle et al.; 2002; 
Hoekstra et al. 2002; Kaiser et al. 2000; Maoli et al. 2001; 
Refregier et al 2002; Rhodes et al. 2001; Van Waerbeke et al. 
2000, 2001, 2002; Wittman et al. 2000). The shape of the cos- 
mic shear signal as a function of angular scale remarkably fol- 
lows theoretical expectations, which successfully confirms the 
gravitational instability scenario, even on small scales where 
non-linear structures dominate the lensing signal. 

Further evidence for lensing is provided by the gravita- 
tional magnification bias. In addition to distortion, distant ob- 
jets are magnified or demagnified, depending on whether the 
matter along their lines-of-sight is over- or underdense com- 
pared to the mean. Magnified sources are preferentially in- 
cluded into flux-limited samples, thus sources behind matter 
overdensities are somewhat over-represented. Since galaxies 
are biased with respect to the dark-matter distribution, it is ex- 
pected that this effect induces cross-correlations between dis- 
tant sources and foreground galaxies. The existence of sig- 



nificant cross-correlations between distant quasars and fore- 
ground galaxies on angular scales of several arc minutes has 
indeed been firmly established (see Bartelmann & Schneider 
2001 for a review) and motivated further theoretical develop- 
ment in order to predict how the magnification bias depends on 
cosmological models. Following earlier work by Bartelmann 
(1995) and Dolag & Bartelmann (1997), Menard & Bartelmann 
(2002) demonstrated the high sensitivity of angular quasar- 
galaxy cross-correlation function to several cosmological pa- 
rameters, namely the matter density parameter, Hq, the nor- 
malisation and shape of the dark-matter power spectrum, as 
and T, and the bias parameter of the galaxies, b. Hence, mag- 
nification bias of quasars is equally efficient as cosmic shear in 
constraining the geometry and the dark matter power spectrum 
of the Universe. However, as with cosmic shear, the informa- 
tion provided by the quasar-galaxy correlation function alone 
is insufficient for independently constraining all these parame- 
ters. 

Following similar motivations as Bernardeau, van 
Waerbeke & Mellier (1997) and Jain & Seljak (1997), we 
decided to explore how deviations from Gaussian statistics 
produced by non-linear growth of structures could modify the 
galaxy-quasar cross-correlation signal and eventually break 
some degeneracies between comological parameters. The eas- 
iest approach is to focus on the additional information that can 
be extracted from higher-order correlations between quasars 
and galaxies which are most sensitive to non-Gaussianity, 
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namely the correlation between distant quasars and foreground 
galaxy pairs. As for the skewness of the convergence field, 
we can expect that some parameter dependencies disappear by 
normalising the three- with two-point correlations. 

The paper is structured as follows. We briefly present 
the formalism of the quasar-galaxy correlation function in 
Sect. ^, assuming the paradigm of gravitational instability of a 
Gaussian random field is valid. In Sect. [}J we then introduce the 
quasar-galaxy-galaxy correlator and predict some useful obser- 
vational signatures. Section ^| deals with density statistics and 
the numerical evaluation of the triple correlator. We then de- 
fine a skewness parameter in Sect. |5] and demonstrate how it 
can be used for directly measuring Slo- Similarly, we show in 
Sect. H how several properties of the galaxy bias can be con- 
strained. We finally estimate the signal-to-noise ratio of the cor- 
responding observation, and specialise it for the Sloan Digital 
Sky Survey in Sect. ||. 

2. Magnification-induced correlation functions 

Basic statistical properties of the magnification due to gravi- 
tational lensing by large-scale structure have been studied in 
earlier papers. To lowest order in the relevant quantities, mag- 
nification is proportional to the lensing convergence, which 
has identical statistical properties to the lensing-induced dis- 
tortions, i.e. the cosmic shear. Those were investigated in many 
studies, starting with the pioneering papers of Gunn (1967) and 
Blandford et al. (1992). 

The two-point correlation function caused by gravitational 
magnification between background quasars and foreground 
galaxies was first introduced by Bartelmann (1995) and gen- 
eralised by Dolag & Bartelmann (1997). We refer the reader 
to these papers for detail and only briefly recall the formalism 
and approximations leading to the two-point correlation func- 
tion. The notation and the definitions used for the three-point 
correlation function are presented in the next section. 

2. 1. The two-point correlation function w(0) 

The angular two-point correlation function between quasars 
and galaxies is defined by 

{[n Q (cj)) - n Q ] [n G (8 + 0_) - n G }) 



(1) 



where uq.g are the number densities of quasars and galaxies on 
the sky. The bar denotes the average on the sky, and the angular 
brackets denote averaging over positions <fi and the directions 
of 0, assuming isotropy. 

Lensing magnification increases the flux received from 
sources behind matter overdensities, but also stretches the sky 
and thus dilutes the sources, modifying their number density. 
The net effect, an increase or a decrease of the source number 
density, is called the magnification bias. It depends on the num- 
ber of sources gained per solid angle by the flux magnification. 
Let a be the logarithmic slope of the source number counts as 
a function of flux, then the number-density fluctuation is 

n Q (0) - (n Q ) 



where 6fi(9) = /i(0) — 1 is the magnification fluctuation as a 
function of position 0. The magnification is related to conver- 
gence k and (the complex) shear 7 by 

A*= [(1 - k) 2 - MT 1 • (3) 
Thus, to first order in k and 7, the magnification fluctuation is 

6fj, = fi - 1 = 2k . (4) 

Assuming that galaxies are linearly biased with respect to 
the dark-matter distribution, we can write 



n G (9) ~ (no) 



bS(9) 



(5) 



where b is the averaged bias factor of all considered galaxies, 
and S is the projection 



5(0) 



dwps(w) S[f K (w)6, w] 



(6) 



of the density contrast 5 along the line-of-sight, weighted by the 
normalised distance distribution ps(w) of the observed galax- 
ies which are cross-correlated with the quasars. Here, w is the 
comoving radial distance along the line-of-sight, fxiw) is the 
comoving angular diameter distance, and the upper integration 
boundary uijj is the comoving radial distance to the "horizon" 
at z — > 00 (For more sophisticated biasing models, e.g. includ- 
ing stochasticity or non-linearity, see Pen 1998; Dekel & Lahav 
1999; Somerville et al. 2001.) 

We can thus write the two-point correlation function intro- 
duced in Eq. ([[]) as 

w QG (8) = b( f i a ~ 1 (cj>)5(<t> + e)) . (7) 

First-order Taylor expansion of the magnification around unity 
yields 



w QG (0) = 26(a-l)(/<0) 5(0 + 0)) 



(8) 



{net) 



(2) 



Like 5, k is a weighted projection of the density contrast along 
the line-of-sight. Specifically, we can write n as 

k(0) = dwp K (w)6[f K (w)0,w], (9) 
Jo 

where p K (w) is the projector 

, a{w') /a'K) 

The ratio between the angular diameter distances fx is the 
usual effective lensing distance, Wq (w) is the normalised dis- 
tance distribution of the sources, in our case the quasars, and 
a{w) is the cosmological scale factor. 

The correlation function w R $ appearing on the right-hand 
side of Eq. (||) can now be related to the power spectrum of the 
density contrast 5, 

w kS (9) = <k(0) 6(4 + 6)) 

= J dwp K (w) J dw' ps(w') 

x (S[f K (w)9,w}S[f K (w')(cj ) + e) 1 w'}) , (11) 
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where the integrations over w and w' range from to u>h- 

We can now use Limber's equation for the statistics of pro- 
jected homogeneous and isotropic random fields. Inserting the 
Fourier transform of the density contrast, and introducing its 
power spectrum Ps(k), we find 



w K s{9) 



dw 7J7 — \ — 

sds 



2ir 



Ps 



fi<(w) : 



>w) J (s9) 



(12) 



3. The three-point correlation function 

3. 1. Formalism 

Let us now extend the formalism introduced in the previous 
section to define higher-order statistical quantities. There are 
two possibilities for defining a three-point correlator between 
quasars and galaxies, either through correlations between sin- 
gle quasars and galaxy pairs, zqgg> or between quasar pairs 
and single galaxies, zqq G . 

By definition, three-point correlations vanish for Gaussian 
random fields. They are created in the course of the non- 
linear evolution of the underlying density field, hence they 
are expected to appear preferentially on small angular scales. 
Since the mean separation between quasars is in general much 
larger than between galaxies, correlations between quasars and 
galaxy pairs should be much easier to measure. We will there- 
fore focus on the triple correlator zqgg only. 

By definition, and using the formalism introduced in the 
previous section, we have 

2qgg(0i, 2 ) = (*M a-1 (^) + 0i) + 2 )} , (13) 

where the average extends over all positions <f>. 

Assuming as before a linear biasing relation between 
the galaxies and the density fluctuations, and expanding the 
Eq. (|l3]) to first order in k and 7 leads to 

zqgg(0i,0 2 ) = 2b 2 (a-l)z(9 1 ,9 2 ) , (14) 
where z(Ox, 6 2 ) is the three-point correlation function 

z(e 1) e 2 ) = {n{ct>)~5{ct> + e 1 )~5{ct> + e 2 )) 

= J dwip K {wi) J dw 2 pg(w 2 ) J dw 3 p s (w 3 ) 
x (S[f K (wi)(l>, wi] S[f K (w 2 )((f> + 0i),w 2 ] 



x 5[f K (w 3 )((l> + d 2 ),W3]) 



(15) 



Again, the integrations over 11)1,2,3 extend from to u>h- 
Stochastic or non-linear biasing schemes will be discussed in 
Sect. |. 

The three-point correlation function Zq GG (8i,8 2 ) is re- 
lated to the excess probability with respect to a random dis- 
tribution for finding triangle configurations, defined by the two 
angular separation vectors 9\ and 9 2 , formed by one quasar 
and two galaxies. As shown by Menard & Bartelmann (2002), 
the lensing-induced quasar-galaxy cross-correlation function 
wqg has a small amplitude, typically on the order of a few 
per cent at angular scales of a few arc minutes. Therefore, in 



order to achieve a significant signal-to-noise ratio for a higher- 
order correlation function such as zqgg> a ver y large number 
of objects will be necessary. Moreover, focusing only on partic- 
ular triangle configurations (for instance equilateral, isosceles 
or any other) will dramatically restrain the number of possible 
measurements in a given survey. For our purpose of measur- 
ing a skewness parameter (see Sect. ||), the detailed angular 
dependence of a given triangle configuration is not of imme- 
diate relevance. Thus, it is observationally preferable to focus 
on an angular average of zq GG (9i,9 2 ) over all suitable tri- 
angle configurations inside a given aperture, thus allowing a 
measurement around each quasar. This point will be detailed 
in the next Section, where we will also show how the averaged 
three-point correlation function can be observed. 

3.2. Observational signature 

The statistical meaning of zq g g(0i,0 2 ) is illustrated by the 
expression for the rms fluctuations in the counts of galaxies 
around a given quasar. Following Fry & Peebles (1980), and 
considering projected (rather than three-dimensional) correla- 
tion functions, the variance of the galaxy counts in cells of solid 
angle S at fixed distance from the quasars can be written 

{N 2 )-{N} 2 = {N)+n 2 G [ d 2 9d 2 9' [w GG (9 - 9') 

Js 

+ z QGG (e,6') -w QG (8)w QG (e')} , (16) 

where the angular brackets denote averages over ensem- 
bles of cells of size S. Thus, we have (N) = n G S + 
n G J s d 2 9 wq G (9), where n G is the galaxy number density on 
the sky. 

If the cells were randomly placed rather than at a fixed po- 
sition relative to a quasar, the variance of the galaxy counts 
would be 

(N 2 ) T -{N) 2 = (N) r + nl f d 2 9d 2 9'w GG (9' -9) ,(17) 

Js 

with (A^) r = n G S. 

The normalised extra variance A(9) of galaxy counts in 
cells of area S near quasars is directly observable and can be 
expressed in terms of the integrated correlation functions intro- 
duced above, 



A(9) 



(N 2 ) - (AO 2 (A 2 ) r - (N)i 



z QGG (9)-w 2 QG (9) + 



(AO? 

WQ G (9) 



(N)r 



(18) 



where wq G and zqgg are the cell-averaged correlation func- 
tions 



^qg(^) 
and 



Jd 2 9'U e (9')w QG (9') 
J d 2 9'U e (9') 



(19) 



- (f) s Jd 2 9 1 U e (9 1 ) Jd 2 9 2 U e (9 2 )z QGG (9 u 9 2 ) 

Zqgg{6) fvoMojfvoMOi) - (20) 
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Here, the function Uq(9') is a top-hat filter of radius 9. Other 
filters can be used if one wishes to be more sensitive to spe- 
cific wavelength ranges of the dark-matter power spectrum (see 
Schneider et al. 1998) but we will not investigate this point in 
the present paper. 

In summary, the relevant third-order quantity for our pur- 
pose is the cell average zqqg(9) of the triple correlator 
zqgg(<?i, #2) between quasars and galaxy pairs. We will now 
estimate its properties and then show how second and third- 
order quasar-galaxy correlations can be combined to measure 
f^o and several properties of the galaxy bias. 



Note that we have expanded the lensing-induced magnifica- 
tion to first order in Eq. (|15| ). The next-order term has a contri- 
bution proportional to K^which, when introduced in Eq. (p2|), 
is of order (5W) 4 and thus formally of the same order as the 
other terms in Eq. (^2|). However, this additional term differs 
by the weight factor /k(w — w')/ from the other terms 
since it contains the lensing efficiency. As a result, this ad- 
ditional factor will be one order of magnitude smaller, even 
though it is of the same order in the perturbation series for 8 
(see Van Waerbeke et al. 2001 for more detail). We can there- 
fore neglect its contribution. 



4. Evaluation of the excess galaxy variance A 

4. 1. Statistics of the density field 

Under the common assumption that the initial density fluctua- 
tions were Gaussian and that cosmic structure grew by gravita- 
tional instability, the three-point correlation function is intrin- 
sically a second-order quantity, and should be detectable only 
where non-linearities arise in the density field. 
We expand the density field to second order as 



(21) 



where 8^ is of order (8^) 2 and represents departures from 
Gaussian behaviour. Note that since two of the three 8 factors 
in Eq.dl^) are derived from the fluctuations in the galaxy num- 
ber density, this expansion is only valid if galaxies are linearly 
biased relative to the dark-matter fluctuations. For now we will 
assume this linearity and show in Sect. || how cosmic magnifi- 
cation can be used to measure the matter density parameter. 
(This property of linearity has been measured on scales be- 
tween 5 and 30/i _1 Mpc with the 2dF Galaxy Redshift Survey 
(Verde et al. 2002)). In next section we will also show that 
quasar-galaxy correlations can also test this property and probe 
the angular range where this simple relation between dark mat- 
ter and galaxies breaks down. 

Assuming linear biasing and expanding 8 to second order, 
we obtain for (818283) 



(818283) 



(4 1) 4 1) 4 1) ) + (s\ 



(1) 4 X) 4 2; 



+ cyclic terms (231,312) 



(22) 



The first term in Eq. (£2|) vanishes because the density fluctua- 
tion field is Gaussian to first order, hence the third moment of 
8^ is zero. Thus, the leading term in Eq. (p2h is of the order 
of^^f). 

In second-order perturbation theory, the Fourier decompo- 
sition of the density-fluctuation field is given by 

8 (2 Hk) = J d 3 kid 3 k 2 S^(ki)8^(k 2 ) 

x 8 B (k! +k 2 -fe)F(fe 1 ,fe 2 ) , (23) 

with 

w(h h \ 5 l fcr-fc 2 fki k 2 \ 2 (fc 1 -fc 2 ) 2 
F(ki,k 2 ) = - + ~ 2l ^ [y2 + Ti) + 7^kT~ ' (24) 



4.2. Bispectrum and Non-Linear evolution 

The ensemble average in Eq. (|22|) is related to the bispectrum 
in Fourier space. By definition, 



(8(ki)8(k 2 )8(k 3 )) = (2nfB(ki,k 2 ,k 3 ) 
x (5 D (fci + k 2 + k 3 ) ■ 



(25) 



Inserting Eq. ( j24| ) into Eq. ( p5| ) leads to an expression of the 
bispectrum in terms of the second-order kernel F(ki, k 2 ) and 
the dark-matter power spectrum 



B(ki,k2,k 3 ) = 2F(ki,k 2 )P(k 1 )P(k 2 ) 
+ 2F(k 2 ,k 3 )P(k 2 )P(k 3 ) 
+ 2F(ki,k 3 )P(k 1 )P(k 3 ) 



(26) 



Note that Eq. (£5j) implies that the bispectrum is nonzero only 
if the wave vectors (ki, k 2 , k 3 ) form closed triangles. 

For describing the bispectrum on all angular scales, we use 
the fitting formula for the non-linear evolution of the bispec- 
trum derived from numerical CDM models by Scoccimarro & 
Couchman (2001), extending earlier work assuming scale-free 
initial conditions. The kernel F(ki, k 2 ) in Eq. ( j26| ) is then sim- 
ply replaced by an effective kernel F (fei, k 2 ), reading 

F eS (ki,k 2 ) = ja(n,ki)a(n,k2) 



2 fki ■ k 2 



c(n, fci)c(n, k 2 ) 



7 \ kik 2 
The coefficients a, b and c are given by 

l + a 8 -°- 2 (z) [0.7Q 3 (n)] 1/2 (g/4)«+ 3 - 5 



(27) 



a(n, k) = 
b(n, k) = 
c(n, k) = 



1 + ( g /4)™+ 3 - 5 
l + 0.4(n + 3) q n+3 

1 + q n+3.5 

1 + 4.5/ [1.5 +{n + 3) 4 ] (2q) n + 3 
1 + (2q) n + 3 - 5 



(28) 



where n = dlogPii n (fc)/dlogA; is the local slope of 
the dark-matter power spectrum at wave number k, q = 
^/^nl(^), and &NL is the nonlinear wave number defined by 
47rfc§, L Piin(&NL) = 1- To be specific, Pu n {k) is the linear 
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Figure 1. Left panel: The two- and three-point quasar-galaxy correlation functions, averaged within disks of radius 6, are plotted 
as a function of angular scale for a flat Universe with f2rj = 0.3. The dashed curve shows the two-point correlation function w(9) 
and the solid line the three-point function z(9). In dotted lines are plotted the corresponding correlation functions for a linear 
evolution of the perturbation. Near one arc minute, the triple correlator z reaches the ~ 2% level, assuming a linear bias of unity. 
At larger angular scales, the amplitude drops steeply. Right panel: We plot the expected extra-scatter of galaxies around quasars 
for a flat Universe with firj = 0.3 (thick line) and considering a galaxy density of 1 per arcmin 2 . We also show the different 
contributions of Eq. fijl. 



power spectrum at the required redshift. Finally, the function 
Qz{n) is given by 



(4 - 2") 
(1 + 2™+!) ' 



(29) 



which is the so-called saturation value obtained in "hyper- 
extended" perturbation theory (HEPT; Scoccimarro & Frieman 
1999). These expressions imply that at large scales, where the 
coefficients a, b and c aproach unity, the tree-level perturbation 
theory is recovered. On the other hand, at small scales, where 
a 2 — > (7/10)Q 3 (ri.)(7^ ' 4 , b — * and c — * 0, the bispectrum 
becomes hierarchical with an amplitude which approximately 
reproduces HEPT for as « 1. For more detail, see Scoccimarro 
& Couchman (2001). 

4.3. Evaluation of zqgg(#) 

The formalism is now in place for computing the triple corre- 
lator z (91,82)- We start from Eq. (|3]) and replace the density 
contrast S by its Fourier transform. Next, we employ the ap- 
proximation underlying Limber's equation, which asserts that 
the coherence length of the density fluctuation field is much 
smaller than the scales on which both projectors p K and p$ vary 
appreciably. Finally, we insert the expression for the bispec- 
trum described in the previous section and find 



,6 2 ) = / dw p K (w) p$(w) 



d2fc i iki-OtfKM 

(2tt) 2 



d 2 fc 2 



Bs(ki,k 2 , —ki - k 2 , w) e 



ik 2 2 fn{w) 



We saw in Sect. 3 that the suitable quantity for observations 
is not directly 2(61,62), but the corresponding cell-averaged 
triple correlator z(9) defined by Eq. (g0|). The angular integra- 
tion yields 



z(9) = 4 / dwp K (w)p s (w) 



d 2 fc! h[ki6f K (w)] 



d 2 fe 2 



(2^) 2 

2 Bs(ki,k 2 , -fei - k 2 ,w) 



kx6f K (w) 
3i[k 2 8fK(w)} 



(30) 



(27T) 2 i! " "~' k 2 8f K (w) 

The observational signature of cosmic magnification is ex- 
pressed in terms of the normalised excess scatter A(9) of galax- 
ies around quasars; see Eq. (18). This quantity involves both 
the two- and three-point cell-averaged correlation functions 
wqg(S) an d ^qgg(^)- For numerically evaluating these ex- 
pressions, we use the description of the non-linear evolution 
of the CDM power spectrum provided by Peacock & Dodds 
(1996) fitting formula. We normalise the power spectrum such 
that the local abundance of galaxy clusters is reproduced, erg = 
0.52 Q-0-52+o.i3n 0) as determined by Eke et al. (1996). 

We assume for simplicity that all quasars are at the same 
redshift z s = 1.5. More realistic quasar redshift distributions 
do not significantly change the following results as long as the 
foreground galaxies are at comparatively low redshift. We ap- 
proximate the redshift distribution of the galaxies by 



ps(z) dz 



(3z< 



T(3//3) 



exp 



-0 



dz 



(31) 
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with f3 = 1.5 and zq = 0.3. 

The slope of the quasar number counts is fairly well con- 
strained by the most recent quasar catalogues. We use the value 
a = 2 suggested by the first SDSS quasar catalogue (Schneider 
et al. 2001) for quasars brighter than 19th magnitude. Finally, 
we assume 6=1 for simplicity. 

The left panel of Fig. [I] shows zqgg(#) (solid line) and 
wqq(9) (dashed line) as a function of angular scale. These 
quantities were computed for a low-density, spatially fiat uni- 
verse (Slo = 0.3, il\ — 0.7). The dotted lines show the 
expected amplitude of wqg(9) if only linear growth of den- 
sity perturbations is taken into account, and the amplitude of 
•zqgg(#) for quasi-linear evolution. The difference between 
the two regimes changes the amplitude of zqgg by approxi- 
mately two orders of magnitude on small angular scales. On 
large scales, the amplitude of zqgg ($) decreases quickly with 
9 since the density field tends to gaussianity as the smoothing 
scale increases, and thus the triple correlator vanishes. 

Interestingly, the amplitude of zqgg is of the order of one 
per cent on arcminute scales. As for the two-point quasar- 
galaxy correlation function, the amplitude and the shape of 
zqgg(#) are ver y sensitive to cosmological parameters. 

In the right panel of Fig. |l[ we plot the measurable quan- 
tity A(6*) which represents the normalised excess scatter of 
galaxies around quasars; cf. Eq. (ff^). Again, we have used the 
ACDM cosmological model, and we assume a galaxy number 
density of 1 arcmin -2 . 

Evidently, zqqq(9) is the main contribution to A(9) on in- 
termediate and large angular scales. Below a few arcminutes, 
the term wqq(9) / (N) T becomes non-neglegible. This contri- 
bution is due to the shot noise of the galaxies, thus this term 
can be lowered when using galaxies with a higher number den- 
sity. 

5. Determining P. and testing the linearity of the 
bias 

The second- and third-order statistics can be used jointly so 
that several parameter dependencies can cancel. The underly- 
ing physical concept is that second-order statistics quantify the 
Gaussian characteristics of a random process, while third-order 
statistics are non-Gaussian contributions. When used together, 
one can in principle measure their relative strength, thus isolat- 
ing those parameters which are most responsible for deviations 
from Gaussianity. 

The reduced skewness (i.e. the ratio of the third and sec- 
ond moments of a distribution) is a useful practical estimate 
of non-Gaussian features in galaxy catalogs. However, in the 
case of cosmic magnification, it is not possible to define the 
skewness in the same way as for cosmic shear with the con- 
vergence field (Bernardeau et al. 1997), since we are not con- 
sidering the autocorrelation properties of a single field, but the 
cross-correlations between two different fields, namely the dis- 
tributions of foreground galaxies and of the lensing conver- 
gence k. Moreover, the angular-averaged three-point correla- 
tion function zqgg discussed in the previous section is not 
symmetric with respect to permutations between quasar and 
galaxy positions. The three-point cross-correlation of quasars 



with galaxy pairs involves the quasar-galaxy cross-correlation 
as well as the galaxy-galaxy auto-correlation. Since the latter 
does not contribute to the two-point quasar-galaxy correlation 
function, we cannot apply the usual mathematical definition of 
skewness which applies to a unique distribution. Instead, we 
define a pseudo-skewness for our purposes by the ratio 

zqgg(0) 



S' 3 (9) = 



QG 



(0) 



(32) 



Much like the skewness of cosmic shear, S' 3 (6) is insensitive 
in the linear regime to the normalisation of the power spectrum 
as, and still remains weakly dependent on erg even in the non- 
linear regime. Moreover, if the galaxy bias can be considered 
linear on certain angular scales, we have 



(33) 



on such scales, where n is close to one half and the ampli- 
tude A{9) can be computed from Eqs. ([jj), (|l3|) and (|l]). 
Interestingly, S' 3 (9) does not depend on the linear galaxy bias 
factor. We plot the pseudo-skewness in Fig. [|for different flat 
and cluster-normalised CDM model universes with density pa- 
rameters ilo = 0.3, 0.5 and 1. The solid line shows the result 
obtained with the non-linear prescription of the bispectrum, 
and the dashed line shows the result of using quasi-linear the- 
ory. 




9 (arcmin) 

Figure 2. The skewness S' 3 (6) is shown as a function of an- 
gle 9 for three different cosmological models. The solid line 
is the prediction based on the non-linear description of the 
power spectrum and the bispectrum. The dashed line shows the 
perturbation-theory calculation. 



The non-linear growth of density perturbations introduces 
some dependency of the pseudo-skewness S' 3 (9) on the nor- 
malisation ag and the shape parameter T of the dark-matter 
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0.6 0.8 1.0 1.2 1.4 0.14 0.16 0.18 0.20 0.22 

Figure 3. The dependence of the skewness S' 3 on erg and T is shown for two different cosmologies (£lo = 1 and 0.3) and three 
different angular scales, 1', 5' and 10'. A value of T = 0.21 is used when as varies and we fix as — 0.9 when T varies. The 
dependence of S' 3 on as goes down to 10% at small scales and is at the one percent level when T changes. 



power spectrum because wqg and zqgg depend on the shape 
of the power spectrum in different way; cf. Eqs. (|2j) and (|j"5|). 
We plot these dependences in Fig. |5} 

The figure shows that S' 3 (9) depends relatively weakly on 
as and is almost insensitive to T. Varying as from 0.5 to 1.4 
changes the amplitude of S' 3 (6) by 30% near 10 and 5 arcmin- 
utes and this change in the amplitude reduces to 10% near one 
arc minute. Therefore we see that the scales of a few arcmin- 
utes will be favored for estimating £lo- The effect of changing 
r from 0.14 to 0.23 is even weaker and does not depend on 
the angular scale as well as the value of Slo- Considering the 
above mentionned range the S' 3 dependence on T is at the one 
percent level. Moreover, if additional measurements allow to 
reduce the possible ranges for as and T, the dependence of S3 
on these parameters further reduces. Thus, we can consider the 
pseudo-skewness S'3 (6>) to be weakly sensitive to normalisation 
of the dark-matter power spectrum, and almost insensitive to its 
shape parameter. Likewise, the cosmological constant 17a has 
a negligible effect. 

On large angular scales, the galaxy bias is expected to be 
linear (Verde et al. 2002), in which case its contribution to S'3 
cancels out. Finally, the only effective parameter we are left 
with is £lo, which means that the pseudo-skewness S' 3 (8) can 
effectively constrain the matter-density parameter. 

The pseudo-skewness S' 3 (9) can also test the linearity and 
stochasticity of the bias parameter which is of interest at 
smaller scales. For a given cosmology and assuming a CDM 
power spectrum, any departure from the angular variation of 
S' 3 (9) given in Fig. ^ can be interpreted as resulting from non- 
linearity and/or stochasticity of the bias. Indeed, our previous 
calculations deriving the two- and three- point correlators used 
only two hypotheses, namely the linearity of the bias and the 



assumption that lensing effects occur in the weak regime, i.e. 
<5/i = 2 k + 0(k 2 , |7 2 |). The latter relation is expected to be ac- 
curate on scales larger than a few arc minutes. Below, medium- 
and strong-lensing effects become non-negligible. Moreover, 
as described in Menard et al. (2002), these effects can be quan- 
tified by expanding the Taylor series of the magnification to 
second order, thus allowing the investigation of smaller scales. 
Therefore, the only remaining explanation for a departure of 
the angular variation of S' 3 (9) is a nonlinearity or stochasticity 
of the biasing scheme. 

The pseudo-skewness can thus probe the angular range and 
the corresponding physical scale where the linear relation be- 
tween dark matter and galaxy fluctuations breaks down. 

6. Expected signal-to-noise ratio 

We now estimate the expected signal-to-noise ratio in mea- 
surements of A(6>). The determination of the effective skew- 
ness also requires a measurement of the two-point quasar- 
galaxy correlation function. We refer the reader to Menard & 
Bartelmann (2002) for a detailed study of the signal-to-noise 
ratio expected for wqg(#)- 

For the measurement of zqqq(9) we ar e n °t aiming at a de- 
tailed noise calculation, but rather an approximate estimation of 
the main source of error, i.e. the finite sampling error caused by 
the limited number of available quasars. As shown in Fig. ^, the 
dominant contribution to the excess scatter A(9) introduced in 
Eq. (^8]) is zqgg(#)i except on small scales. We therefore use 
the simplifying assumption that A ~ zqgg m the following. 

Since the excess scatter of galaxies around quasars defined 
in Eq. (|l8]) is a counts-in-cells estimator, its measurement accu- 
racy will be limited by the finite size of the available sample, by 
boundary and edge effects, and by the effects of discrete sam- 
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pling. In practice, measurements of zqgg($) will t> e restricted 
to angular scales much smaller than the size of the survey, thus 
the errors contributed by boundary effects can be considered 
negligibly small compared to the finite sampling of the galax- 
ies, which causes the main limitation. 

In order to estimate this noise, we assume the quasars to be 
at random positions in the sky, i.e. uncorrelated with the galaxy 
positions. In fact, physical correlations are excluded given the 
required large separation between the two populations, and 
on the other hand the cross-correlations between quasars and 
galaxies induced by lensing are so weak that the correspond- 
ing change in the galaxy distribution is entirely negligible in 
the total error budget. We then assume a realistic distribution 
of galaxies (see Appendix A) and find the standard deviation 
of the normalised scatter of the galaxy counts around Nqso 
quasar positions to be 



(N 2 ) - (N)' 



— (- 

N Q so VAT 3 

8Aw 2 



N 



2w 



2 + 3w 
N 2 

-28.2w) 3 



Aw 

~w 



(34) 



where N is the average number of galaxies in cells of a ra- 
dius 8. This result is derived in Appendix A, using the formal- 
ism developed by Szapudi & Colombi, and considering galaxy 
properties measured by SDSS. 




6 (arcmin) 

Figure 4. Number of quasars required for a 3-a detection of the 
excess scatter A of the galaxy counts near quasars, as a func- 
tion of angular scale 6, and according to our simplified noise 
estimation. Three different values for the galaxy number den- 
sity are assumed, namely 1, 2, and 3arcmin~ 2 . The vertical 
line indicates the mean angular separation of 20' between dis- 
tant quasars in the SDSS. 



The signal we are interested in is an excess scatter in galaxy 
counts. Therefore, the noise of the total measurement is twice 
the value of a introduced in Eq. ([34|). It is possible to re- 
duce this noise by measuring the scatter of the galaxy counts 
at a larger number of random locations since there are typi- 
cally many more available galaxies than quasars, therefore the 
dominant contribution comes from the measurements around 
quasars only. Finally, we can write the number of required 
quasars for achieving a detection at an angular scale : 



Nqso — 



[zqgg(0) ere ] 



(35) 



where v is the desired signal-to-noise ratio and erg is the error 
computed from Eq. 34 We recall that the quasar number given 
by Eq. ( ^Bj ) refers to quasars with redshifts higher than those of 
the galaxy population. 

Several weighting schemes can be used to maximise the 
signal-to-noise ratio of the detection. Menard & Bartelmann 
(2002) showed how to optimally weight the contribution of 
each quasar with respect to its magnitude. Weighting the galax- 
ies with respect to their redshift can also increase the signal- 
to-noise ratio of the detection. Indeed the triple correlator 
zqgg(0i , #2) is related to the excess of triangle configurations 
in which two galaxies trace a high density region of the dark 
matter field, i.e. configurations in which the two galaxies are 
close in angle and redshift. Projection effects mimicking galaxy 
pairs will contribute noise to the final measurement. Therefore, 
the width of the galaxy redshift distribution inside a given cell 
around a quasar can be used to additionally weight the final 
measurement. Numerical simulations will be needed to quan- 
tify the change of the final signal-to-noise ratio, as well as the 
effects of galaxy clustering and cosmic variance which were 
not taken into account in our noise estimation. 

As an application, we now investigate the feasibility of 
measuring A(9) with the data of the Sloan Digital Sky Survey 
(SDSS ; York et al. 2000). Within this project, the sky has 
already been imaged for two years, and the survey will be 
completed in 2005, reaching a sky coverage of ~ 10 000 
square degrees. Depending on the limiting magnitude of the 
selected sample, SDSS can achieve a galaxy density of ?ig ~ 
1 arcmin -2 for galaxies observed down to r' = 21, or uq 
3 arcmin -2 down to r' = 22, but requiring extensive careful 
masking of regions with poor seeing within the survey, and a 
careful star-galaxy separation (Scranton et al. 2001). 

In Fig. |J, we present for different values of the galaxy num- 
ber density the number of quasars required in order to achieve 
a 3-er detection of the excess scatter A(0). The figure shows 
that the measurement becomes more easily reachable at inter- 
mediate angular scales. Note that in our estimation we have 
assumed isolated quasars and Eq. [35] ceases to be valid at angu- 
lar scales reaching the average angular separation of quasars. In 
the case of SDSS, this angular scale is close to 20' for the sam- 
ple of spectroscopic quasars satisfying z > 1. On larger angu- 
lar scales, the real errors will increase due to correlated galaxy 
counts, and the number of quasars estimated from Eq. ( |35] ) will 
no longer be reliable. The SDSS will observe 10 5 spectroscopic 
and 5 x 10 5 to 10 6 photometric quasars. Thus, we conclude 
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from Fig. [| that SDSS will exceed the number of quasars re- 
quired for measuring three-point correlations between quasars 
and galaxies, thus allowing a new direct measurement of the 
matter density f2 - 

7. Conclusion 

Via the magnification bias, gravitational lensing by large-scale 
structures gives rise to angular cross-correlations between dis- 
tant sources and foreground galaxies although the two pop- 
ulations are physically uncorrelated. Depending on whether 
the matter along the lines-of-sight towards these background 
sources is over- or underdense with respect to the mean, and 
depending on the value of the slope a of the cumulative num- 
ber counts of the sources, magnification effects can cause an 
excess or a deficit of distant sources near foreground galaxies. 

These lensing-induced correlations carry information on 
the projected dark-matter distribution along the lines-of-sight 
and can thus provide constraints on cosmology. Considering 
distant quasars and foreground galaxies, Menard & Bartelmann 
(2002) quantified these constraints and showed that given the 
large number of parameters involved (the matter density pa- 
rameter, f2o> the normalisation and shape of the dark-matter 
power spectrum, as and T, respectively, and the bias parame- 
ter of the galaxies, 6), the information provided by the quasar- 
galaxy correlation function alone is insufficient for indepen- 
dently constraining all of these parameters. 

In this paper, we have investigated what additional informa- 
tion can be expected from higher-order statistics. Such statistics 
have similar weightings of the power spectrum along the line- 
of-sight and can measure non-Gaussianities of the density field 
due to the non-linear growth of structures. We have specifically 
considered correlations beween distant quasars and foreground 
galaxy pairs and showed that this three-point correlator can be 
related to the excess scatter of galaxies around quasars with 
respect to random positions, which is a straightforwardly mea- 
surable quantity. Using the assumptions that 

- galaxies are linearly biased with respect to the underlying 
dark matter, 

- the dark-matter distribution is described by a CDM power 
spectrum, 

- and the lensing magnifications can be approximated by /i = 
1 + 2k in the weak lensing regime, 

we have computed the amplitude of the expected excess scatter 
of galaxies in the vicinity of quasars. For a flat cosmology with 
f2rj = 0.3 and as = 0.93, we find the amplitude of this effect to 
reach the per cent level at angular scales near one arc minute. 

We further showed that combining second- and third-order 
statistics allows a pseudo-skewness parameter S 3 to be defined 
which turns out to be weakly sensitive to the normalisation and 
the shape of the power spectrum. Moreover, if the linear biasing 
scheme is valid, this parameter is only sensitive to the matter 
density flo; the dependences on the other parameters (erg, T, A 
and 6) are weak or cancel out completely. Thus the skewness 
S 3 provides a direct and independent measurement of firj- 

We computed the expected angular variation of S'3 and 
showed that for a given cosmology and assuming a CDM power 
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spectrum any departure from the predicted angular shape must 
be due to a non-linear and/or stochastic behaviour of the galaxy 
bias, on scales larger than a few arc minutes. 

Finally, we estimated the signal-to-noise ratio of the ex- 
pected excess scatter of galaxies near quasars, which is the 
main source of noise in S'3. We derived the noise coming from 
the finite sampling of galaxies having realistic distributions. 
Applying our result to the Sloan Digital Sky Survey, we find 
that S 3 should be measurable on large scales from that survey 
with about twenty thousand distant and bright quasars, allow- 
ing thus a direct and independent measurement of f2o- 

On smaller angular scales less quasars are required and the 
parameter S 3 can probe the angular range and the correspond- 
ing physical scales where the linear relation between dark mat- 
ter and galaxy fluctuations may break down. 
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Appendix A: Noise estimation 

We use the formalism developed by Szapudi & Colombi (1996) 
for calculating errors due to finite sampling, and we refer 
the reader to this paper as well as to the review on "Large- 
Scale Structure of the Universe and Cosmological Perturbation 
Theory" by Bernardeau et al. (2001). We first define some use- 
ful quantities and then derive the error on our estimator due to 
Poissonian noise. 
We have the estimator 

Let Pjy denote the probability of finding N galaxies in a cell of 
a given size. The factorial moments are defined as 

F k = ((N) k )=J2(N)kP N , (A.2) 

where (N) k = N(N - l)...(N - k + 1) is the fc-th falling 
factorial of N, and the ensemble average can be evaluated using 
the known probability distribution Pjy. We can write the error 
on our estimator in terms of factorial moments, 

where the operator (. . .)c averages over all possible ways of 
throwing C cells into the survey volume. In order to compute 
the difference, we now perturb F\ and F 2 as F\ = F\(l + £=r*) 

and F 2 = F 2 (l + ^), with (£Fi) = (SF 2 ) = 0. 
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This yields 



E z 



4F 2 



4£f 



(SFl) 



2F[ 



*0.(5F 1 6F,) + *£i.(8F*). 



Ft 



(SFi 5F 2 ) 



(A.4) 



For evaluating this expression, we need to compute the vari- 
ances and the co variance of F\ and F2. We will do this using 
the generating f uncti on of the probability distribution Pm in- 
troduced in Eq. (A. 2), 



p(x)=Y i p i 



NX 



.,N 



(A.5) 



The nice property of this quantity is that the factorial moments 
are obtained by a Taylor expansion of P(x) around x = 1. The 
cosmic covariance on factorial moments can then be written as 



x—y—l ' 



(A.6) 



and for estimating the finite-sampling error, the error generat- 
ing function £ is 



£(x,y) 



P{xy) - P(x)P(y) 
C 



(A.7) 



(Szapudi & Colombi 1996), where C is the number of cells 
used, i.e. the number of quasars or random positions involved 
in the measurement. Thus, we can compute the variances and 
covariances needed in Eq. ( A.4 ), 



1 

c 
1 

c 



N(N - 1) P N + N P N - (J2 N Pn) 2 



L N 



N 



(F 2 +F 1 -Ff] 



(A.8) 



For Poissonian noise, we have the simple relation F k = N k , 
where N is the average number of galaxies in a cell. Thus, the 
expression reduces to 



In the same way, we obtain 



C dx 2 dy 2 
1 

C 



Y^N 2 (N- l) 2 P N - (£2 N(N - 1) P N ) 



N 



N 



i(F 4 + 4F 3 + 2F 2 -F 2 2 ) 
47V 3 + 2iV 2 



1 

C 



N 2 (N - 1) Pjv - J] N(N - 1) P N 



N 



N 



N 
1 



c_ 

2N 2 



(F3 + 2 F2 — F\ F2) 



(A.ll) 



Inserting those terms into Eq. (A.4), we find after further 
simplifications 



E' 



( (N 2 )~(Nf 
\ (N) 2 



1 



C N 2 



CN 3 



(A.12) 



which is the error due to the finite quasar number. 

In order to check this result, we have performed a numer- 
ical simulation in which we measured our estimator on two- 
dimensional Poisson distributions of particules. We then com- 
puted its standard deviation for several cell numbers and partic- 
ular point densities and find the numerical results in full agree- 
ment with the previous expression. 

The clustering of galaxies can also be taken into account 
in this calculation. To do so, we express the factorial moments 
of the galaxy number counts in terms of n-points cell-averaged 
correlation functions : 



F k = N k ((1 + 5) fe ) . 



(A.13) 



Using the fact that (<5 2 ) = (S 2 ) c = w{9), (S 3 ) = {6 3 ) c = z{6) 
and (c5 4 ) = (<5 4 ) c + 3 (<5 2 ) 2 = w 4 (6>) + 3 w{9) 2 where w, z and 
W4 are the 2-, 3- and 4-point cell-averaged correlation functions 
of the galaxies and the subsript c denotes the connected contri- 
bution, we can use the expression of the factorial mome nts i n 
terms of correlation functions in the error expression (Eq. A.4 ). 
One finds : 



E 



2 f (N 2 )-(N) 2 
V (N) 2 



1 



2 + 3w 4w + 2z 



CN 3 CN 2 



CN 



W4 - 4zw + 2w + 4tD J . (A.14) 



Under the hierarchical assumption, higher-order moments can 
be computed by using the skewness and kurtosis parameters of 
(A. 9) tne g a j aX y distribution : z = S3 w 2 and W4 = S4 u> 3 . The pre- 
vious expression can therefore be expressed in terms of w(9) 
only. Then, each term can be evaluated by considering a power- 
law correlation function : w(9) = 9q ( 1 ^ and using 
SDSS measurements of galaxies with 20 < r < 21 in order to 
estimate each parameters. We have : 9$ = 1CP 2 ' 3 and 7 = 0.7 
(Connolly et al. 2002) and for the higher-order moments of the 
same galaxies : S3 = 4.2 ± 0.4 and S4 = 41 ± 10 (Szapudi et 
al. 2002). 



and 

{SF 1 SF 2 } 



C 

1 d d 2 
C dx dy 
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